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Safe  and  reliable  operation  of  a  Li-ion  battery  requires  control  and  often  management  of  the  thermal 
envelope.  In  this  context,  a  two-dimensional,  transient  mathematical  model  comprising  conservation  of 
charges,  species,  and  energy  together  with  electroneutrality,  constitutive  relations  and  relevant  initial  and 
boundary  conditions  for  a  spiral-wound  cylindrical  Li-ion  battery  is  derived  and  solved  numerically  for 
passive  thermal  management  with  and  without  a  phase  change  material  (PCM)  at  various  galvanostatic 
discharge  rates.  Two-way  coupling  of  the  electrochemical  and  thermal  equations  of  change  is  attained 
through  heat  generation  terms  and  temperature-dependent  physical  properties.  Within  this  framework, 
the  electrochemical  and  thermal  behavior  is  discussed  in  terms  of  edge  effects  arising  from  the  design 
of  the  spiral-wound  structure  and  variations  in  heat  generation  in  the  functional  layers.  In  addition, 
the  cell  performance  with  passive  thermal  management  through  a  PCM  is  shown  to  lower  the  overall 
temperature  of  the  cell  at  discharge  rates  around  5  C-rates,  provided  the  PCM  layer  is  thick  enough  to 
provide  cooling  during  the  entire  discharge.  The  model  can  be  employed  for  wide-ranging  parameter 
studies  as  well  as  multi-objective  optimization  of  not  only  design  parameters  pertaining  to  the  spirals 
but  also,  for  example,  for  design  of  the  thickness  of  the  PCM  layer. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Recent  years  have  seen  several  recalls  of  commercial  Li-ion  bat¬ 
teries  due  to  overheating.  The  heat  originates  mainly  from  the 
electrochemical  reactions  that  occur  during  charge/discharge  of 
the  battery  as  well  as  Joule  heating  -  if  this  heat  is  not  dissi¬ 
pated  properly,  it  can  lead  to  overheating  of  the  battery,  and, 
in  a  worst-case  scenario,  thermal  runaway.  The  latter  typically 
occurs  in  conjunction  with  ‘abuse’  conditions,  high  power  draw, 
and/or  manufacturing  defects.  Generally,  some  form  of  protection 
is  thus  needed  to  prevent  overheating,  either  by  means  of  elec¬ 
tronic  control  circuits  that  monitor  charge/discharge  rates  and/or 
temperature,  or  by  means  of  a  properly  designed  thermal  man¬ 
agement  system.  Furthermore,  for  the  purpose  of  operating  the 
battery  within  safe  operating  temperature  limits,  either  one  of 
two  following  thermal  management  strategies  are  typically  intro¬ 
duced:  active  systems  with  air/liquid  cooling  [1  ]  or  passive  systems 
with,  for  example,  phase-change  materials  (PCMs)  [  1  -6] ;  see  Band- 
hauer  et  al.  [7]  for  a  detailed  review  of  these.  The  main  advantages 
of  passive  cooling  with  PCM  viz-a-viz  active  cooling  are  a  sim¬ 
plified  design,  absence  of  parasitic  power  consumption,  smaller 
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temperature  gradients  with  air  as  coolant  under  normal  and 
stressed  operating  conditions  [1],  and  that  propagation  of  energy 
from  cell  to  cell  arising  from  thermal  runaway  inside  a  battery  pack 
can  be  reduced  [5]. 

Various  mechanistic  mathematical  models  have  been  devel¬ 
oped  in  order  to  predict  the  transient  electrochemical  and  thermal 
behavior  of  a  Li-ion  cell  in  a  rectilinear  geometry  [8-17]  and  a 
spiral-wound  geometry  [18-25].  In  essence,  these  models  typically 
consider  the  transient  equations  of  change  for  species,  charge  and 
energy  together  with  relevant  boundary  conditions  and  constitu¬ 
tive  relations.  The  level  of  detail  and  resolution  depends  on  the 
treatment  of  electrochemical  and  thermal  phenomena:  the  most 
detailed  models  [11,13-15]  solve  the  governing  equation  in  the 
form  of  partial  differential  equations  (PDEs),  thus  resolving  the 
local  transport  phenomena,  electrochemistry  and  heat  generation. 
Loss  of  detail  is  often  incurred  when  one  of  the  governing  equa¬ 
tions  is  simplified  to  an  ordinary  differential  equation  (ODE)  or  a 
phenomenological  expression,  thus  only  truly  accounting  for  the 
global  behavior  of  the  battery  cell  and,  to  a  lesser  extent,  the  local 
behavior;  for  example,  when  only  the  equation  of  change  for  energy 
is  solved  locally  together  with  an  expression  for  heat  generation 
from  some  form  of  approximative  relation  [17-19,21,22,24]  of  the 
electrochemical  reactions  or  from  experimental  measurements 
[  1 ,5  ] ;  or  when  the  equations  of  change  describing  the  electrochem¬ 
ical  phenomena  on  a  local  level  are  coupled  with  some  form  of 
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Nomenclature 

a*  constants  in  electrolyte  conductivity  expression 

A  specific  interfacial  area  per  unit  volume  (m-1 ) 

Q  electrode  capacity  (Ah  nrr 2 ) 

Cth  theoretical  capacity  of  electrode  material  (mAh  g-1 ) 

C[  electrolyte  concentration  ( mol  m-3 ) 

cf°  initial  electrolyte  concentration  (mol  nrr3 ) 

Cp  effective  specific  heat  capacity  (J  (kg  K)-1 ) 
c°  initial  concentration  of  lithium  in  the  active  material 

(molm-3) 

cs  concentration  of  lithium  in  the  active  material 

(molnrr3) 

cf vg  average  concentration  of  lithium  in  the  active  mate¬ 

rial  (mol  nr3) 

csmax  maximum  concentration  of  lithium  in  the  active 
material  (mol  nr3) 

cssurf  surface  concentration  of  lithium  in  the  active  mate¬ 
rial  in  the  electrodes  (mol nr3) 

D/  diffusion  coefficient  of  electrolyte  (m2  s-1 ) 

Ds  diffusion  coefficient  of  lithium  in  the  active  material 

(m2  s-1) 

Ea  q  activation  energy  for  a  variable  ©  (kj  mol-1 ) 

F  Faraday’s  constant  (96,487  C  mol-1 ) 

h  height  of  the  battery  (m) 

f)  heat  transfer  coefficient  (W  nr2  K-1 ) 

iapp  applied  current  density  (A  nr 2 ) 
i'o  exchange  current  density  (A  m-2 ) 

i /  liquid  phase  current  density  (A  irr 2 ) 

is  solid  phase  current  density  (A  nr 2 ) 

if  faradaic  transfer  current  density  (A  m-2 ) 

J  local  charge  transfer  current  per  unit  volume 

(Am-3) 

k  effective  thermal  conductivity  (W (m  K)-1 ) 

k0  reaction  rate  constant 

L  latent  heat  of  PCM  (J  kg-1 ) 

C  length  of  the  spiral  (m) 

ls  diffusion  length  (m) 

Wj  thickness  of  the  layer  i  (m) 

N /  species  (lithium  ion)  flux  (mol  m-2  s-1 ) 

n  normal  vector 

ni,  pi  constants  in  the  entropic  heat  term  for  negative  and 
positive  electrodes 

i%  i<p  constants  in  the  open  circuit  potential  for  negative 
and  positive  electrodes 
Q  volumetric  heat  generation  ( W  nr3 ) 

q  conductive  heat  flux  ( W  nr 2 ) 

R  gas  constant  (J  (mol  K)-1 

Rb  radius  of  battery  (type  1 8650)  (m) 

radius  of  active  material  in  the  electrodes  (m) 
t  time  (s) 

£+  transference  number  of  cation 

T  temperature  (I<) 

Ts  start  temperature  of  phase  change  (I<) 

7/  end  temperature  of  phase  change  (K) 

Urefi  open  circuit  potential  of  the  electrode  i  (V) 
p  number  of  moles  of  Li  in  C6 

t)  number  of  moles  of  Li  in  Mn204 

Greek  letters 

oia  anodic  transfer  coefficient 

ac  cathodic  transfer  coefficient 

£/  volume  fraction  of  electrolyte 

£j  volume  fraction  of  conductive  filler  additive 


£p  volume  fraction  of  polymer  phase 

1 7  overpotential 

0ne  state  of  charge  of  negative  electrode 
Ope  state  of  charge  of  positive  electrode 
v+,  V-  number  of  cations  and  anions  into  which  a  mole  of 
electrolyte  dissociates 

§  emissivity  of  the  outer  can  material 

p  effective  density  (kg  m-3 ) 

a  Stefan-Boltzmann  constant  (W  m-2  K-4) 

( 7/  ionic  conductivity  of  electrolyte  (S  m-1 ) 

ex s  electronic  conductivity  of  solid  matrix  (S  m-1 ) 

0/  liquid  phase  potential  (V) 

0s  solid  phase  potential  (V) 

0°  initial  liquid  phase  potential  (V) 

0°  initial  solid  phase  potential  (V) 

0(T)  placeholder  for  a  temperature  dependent  property 

Subscripts 
amb  ambient 

cc  current  collector 

el  electrolyte 

ne  negative  electrode 

pe  positive  electrode 

ref  reference  value 

sp  separator 

Superscripts 
0  initial  value 

eff  effective  value 


lumped-parameter  model  for  the  thermal  part  on  the  global 
level  [25,26].  The  general  trend  here  is  that  detailed  local  models 
have  been  employed  for  rectilinear  geometries  and  simplified 
counterparts  for  spiral-wound  geometries  -  the  latter  is  typically 
significantly  more  expensive  to  solve  from  the  numerical  point- 
of-view  as  can  be  inferred  from  intricate  geometrical  features  in 
Fig.  1. 

In  view  of  the  lack  of  detailed,  local  resolution  for  modeling  and 
simulation  of  coupled  electrochemistry,  transport  phenomena  and 
heat  generation  in  spiral-wound  geometries,  the  aim  of  this  work 
is  twofold:  First,  to  develop  a  coupled  thermal-electrochemical 
model  for  a  cylindrical  spiral-wound  lithium-ion  battery  with¬ 
out  compromising  local  resolution,  which  can  easily  be  applied 
to  spiral-wound  prismatic  cells  as  well;  and  second,  to  apply  the 
derived  model  to  investigate  the  design  and  operation  of  a  passive 
thermal  management  system  based  on  PCM.  In  short,  the  model 
considers  transient  conservation  of  charges,  species  and  energy; 
it  couples  the  electrochemical  and  thermal  behavior  through  the 
heat  generation  arising  from  reversible,  irreversible  and  ohmic 
heating  as  well  as  through  the  temperature-dependent  transport 
and  electrochemical  parameters.  The  PCM,  in  turn,  is  wrapped 
around  the  battery  and  solved  for  in  terms  of  conservation  of 
energy.  The  results  are  discussed  with  emphasis  on  transient 
behavior  and  temperature  distribution  in  the  various  layers  of 
the  spiral-wound  battery  under  galvanostatic  discharge  at  various 
rates. 

2.  Mathematical  formulation 

We  consider  a  commercially  available  spiral-wound  cylin¬ 
drical  Li-ion  18650  battery  (Fig.  la),  for  which  we  resolve  a 
two-dimensional  cross-section  of  the  battery,  as  illustrated  in 
Fig.  lb,  where  the  functional  layers  -  positive  electrode  (pe), 
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Positive  electrode  (pe) 
Negative  electrode  (ne) 


^Current  collector  (cc) 
^Electrolyte  (el) 


□  Separator  (sp)  |  Outer  can 
|  PCM  Q  Active  material 


Fig.  1.  Schematic  of  (a)  a  1 8650  Li-ion  battery,  (b)  cross-section  of  the  1 8650  battery  showing  the  spiral-wound  jelly  roll,  (c)  cross  section  of  the  battery  with  PCM  surrounding 
it,  (d)  various  functional  layers  in  the  jelly  roll  with  the  roman  numerals  indicating  the  interfaces  of  these  layers  at  the  inner  end  of  the  spiral,  (e)  outer  end  of  the  spiral 
with  the  interfaces  of  various  layers  shown  by  the  roman  numerals,  (f)  modified  computational  domain,  see  numerics,  (g)  agglomerate  structure  in  the  negative  electrode 
(*positive  electrode  also  exhibits  similar  structure)  and,  (h)  diffusion  of  lithium  in  active  material  in  the  electrodes  on  the  microscale. 


negative  electrode  (ne),  current  collector  (cc)  and  separator  (sp) 
-  are  wound  up  in  the  form  of  a  jelly  roll.  We  have  taken  the 
dimensions  of  the  various  functional  layers  from  a  Sony  cell  [27] 
with  the  number  of  wounds  determined  to  be  15  based  on  the 
diameter  of  the  battery  and  the  thickness  of  the  wounds. 

Generally,  the  porous  electrodes  consist  of  active  material, 
conductive  filler  additive,  binder  and  liquid  electrolyte  (el).  In 
this  study,  the  active  material  in  the  negative  electrode  is  LipC6 
and  in  the  positive  electrode  it  is  Lit)Mn204.  The  salt  is  LiPF6 
in  a  nonaqueous  1:2  liquid  mixture  of  ethylene  carbonate  and 
dimethyl  carbonate.  The  electrochemical  reactions  that  occur  at 


the  electrode/electrolyte  interface  during  charge  and  discharge 
are  then 

Discharge 

LipC6  ^  C6  +  ?Li+  +  ?e_ , 

Charge 

Discharge 

LiH_rMn204  +  pLi+  +  pe-  <=*  Lit}Mn204, 

Charge 

where  p  is  the  stoichiometric  coefficient  or  the  number  of  moles  of 
lithium  present  in  the  graphite  structure,  C6,  and  t)  is  the  number 
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Table  1 

Parameters. 


Parameter 

Unit 

cc(-) 

ne 

sp 

pe 

cc(+) 

Reference 

c° 

mol  m3 

- 

2x  103 

- 

[10] 

Cp 

J(kgK)-1 

3.8  xlO2 

7.0  x  102 

7.0  x  102 

7.0  x  102 

8.7  x  102 

[15] 

Cs 

mol  nrr3 

- 

1.5  xlO4 

- 

3.9  x  103 

- 

[10] 

rmax 

mol  ITT3 

- 

2.6  x  104 

- 

2.3  xlO4 

- 

[10] 

D, 

m2  s_1 

- 

7.5  xlO”11 

- 

[10] 

Ds 

m2  s_1 

- 

3.9  x  10-14 

- 

1.0  xlO”13 

- 

[10] 

Ea,D, 

kj  mol-1 

- 

10 

- 

[15] 

Ea,Ds 

kj  mol-1 

- 

4 

- 

20 

- 

[15] 

Ea,at 

kj  mol-1 

- 

20 

- 

[15] 

h 

m 

CTi 

O 

X 

o 

- 

W  ITT2  K-1 

5 

- 

iapp  (1C  rate) 

Am-2 

4.5  x  105 

- 

k 

W(mK)-1 

3.8  xlO2 

0.05  x  102 

0.01  x  102 

0.05  x  102 

NJ 

o 

X 

o 

[15] 

Rb 

m 

9  x  10-3 

- 

m 

m 

- 

12.5  x  10“6 

- 

8.5  xlO”6 

- 

[10] 

Tamb>  Tref 

K 

298.15 

[15] 

Wj 

m 

18  x 10-6 

00 

00 

X 

o 

25  x 10-6 

00 

o 

X 

o 

25  x 10-6 

[27] 

OiayOic 

- 

- 

0.5 

- 

0.5 

- 

[10] 

Sp 

- 

- 

0.14 

- 

0.19 

- 

[10] 

Si 

- 

- 

0.36 

0.72 

0.44 

- 

[10] 

Sf 

- 

- 

0.03 

- 

0.07 

- 

[10] 

q 

- 

- 

0.56 

- 

0.17 

- 

[10] 

p 

kg  m-3 

9.0  x  103 

1.9  x  103 

1.2  xlO3 

4.1  x  103 

2.7  xlO3 

[10,15] 

CTS 

Sm-1 

o 

X 

o 

1  x  102 

- 

3.8 

3.8  xlO7 

[15] 

of  moles  of  lithium  in  the  spinel  structure  of  manganese  dioxide 
Mn204;  Li+  is  the  lithium  ion. 

The  materials  for  the  positive  and  negative  current  collectors 
and  the  outer  can  are  aluminum,  copper  and  stainless  steel,  respec¬ 
tively. 

The  PCM,  which  is  coated  around  the  battery  (Fig.  lc),  is  taken 
to  be  paraffin  wax  impregnated  in  a  graphite  matrix.  Paraffin  wax 
has  a  high  latent  heat  but  a  low  thermal  conductivity,  whence 
the  graphite  matrix  is  provided  to  enhance  the  heat  transfer  rate 
between  the  cell  and  the  ambient  by  conduction  [1,28].  Overall, 
PCMs  have  the  advantage  of  storing  and  releasing  heat  within  a 
narrow  temperature  range  as  latent  heat.  However,  most  of  the 
PCMs  have  a  low  thermal  conductivity  [29]  and  various  attempts 
have  been  carried  out  to  improve  the  conductivity  and  in  turn,  the 
efficiency  of  the  thermal  energy  storage  [28].  In  addition,  factors 
like  mechanical  strength  and  electrical  properties  will  constrain 
the  thickness  and  type  of  the  PCM  that  can  be  employed. 

The  length  scales  ranging  from  the  agglomerate  level  on  the 
order  of  10-7m  to  the  cell  level  on  the  order  of  10_1m  are 
resolved:  the  transport  in  the  cell  on  length  scales  larger  than  that 
of  the  agglomerates  (Fig.  lg)  is  referred  to  as  the  transport  at  the 
macroscale ,  which  includes  mass  transfer  in  the  electrolyte  describ¬ 
ing  the  movement  of  mobile  ionic  species,  material  balances, 
current  flow  and  electroneutrality  based  on  concentrated-solution 
theory  [30],  electronic  charge  conduction  in  the  solid  phase  and 
energy  transfer  in  the  solid/liquid  phases  (Fig.  lb-f);  the  diffusion 
of  ions  in  the  active  material  in  the  electrodes  is  referred  to  as  the 
transport  at  the  microscale ,  which  includes  diffusion  of  lithium  in 
the  active  material  of  the  porous  electrodes  (Fig.  lg). 

The  main  postulates  and  features  of  the  model  are  as  follows: 

1.  Reduction  in  dimensionality.  The  3D  battery  (Fig.  la)  is  reduced 
to  a  2D  cross-section  (lb)  through  the  middle  of  the  battery, 
which  is  justified  by  the  following  arguments:  first,  insulated 
conditions  for  the  energy  flux  are  prescribed  at  the  top  and 
bottom  surface  of  the  battery  (z-direction  in  Fig.  la),  which  is 
a  common  assumption  [19,31,32,18,22,25].  Second,  as  a  first 
approximation,  we  assume  that  the  air  temperature  is  constant 
around  the  battery,  whence  the  temperature  difference,  A Tz, 
inside  the  battery  in  the  axial  direction  is  zero;  i.e.,  ATz  =  0. 
Third,  potential  losses  in  the  axial  direction  in  the  current 


collectors  are  negligible  at  leading  order,  since  the  potential 
drop  in  the  z-direction,  A 0SiZ  ~  iapphl<7s  ~  1 0-4  V  <  Ecen  ~  1 V,  for 
typical  operating  and  material  properties:  iapp  ~  106  Am-2  (~1 
C-rate),  10-3  m,  crs~  107  S  nrr1,  where  iapp  is  the  applied  cur¬ 
rent  density,  h  is  the  height  of  the  battery,  os  is  the  electrical 
conductivity  and  Eceu  is  the  cell  voltage.  Fourth,  the  placement  of 
current  collecting  tabs  on  the  behavior  of  the  cell  is  assumed  to  be 
negligible. 

For  the  resulting  2D  cross-section,  we  prescribe  the  cur¬ 
rent  as  entering  from  the  innermost  boundary  of  the  current 
collector  (II)  and  leaving  at  the  outermost  current-collector 
boundary  (V),  which  is  justified  by  the  potential  drop  in  the 
tangential  direction,  A 0Sft,  around  the  total  wound  of  length, 
£,  which  is  negligible  at  leading  order  compared  to  the  over¬ 
all  cell  voltage;  i.e.,  A (f)s,t~iapp£/<Js  ~  10-2V  «  £ce/|~l  V,  for  a 
typical  length  of  a  wounded  layer,  £~1 0-1  m.  Note  that  this  con¬ 
dition  leads  to  errors  when  the  C-rate  is  around  10  or  higher, 
since  iapp  ~  107  Am-2 (~  10 C-rate),  such  that  A0s,t~iappC/crs  ~ 
10-1  V  <  Eceu~\  V.  Here,  we  limit  the  discussion  to  C-rates  from 
one  to  five. 

2.  Natural  convection.  Natural  convection  in  the  type  of  PCM  that 
we  consider  here  has  been  shown  to  be  negligible  [33].  Inside 
the  battery,  we  postulate  negligible  natural  convection  as  a  first 
approximation  by  letting  gravity  act  in  the  z-direction,  in  which 
there  are  no  temperature/concentration  gradients  due  to  the 
conditions  outlined  in  #1  above;  we  note,  however,  that  one 
could  extend  the  model  to  include  natural  convection  [34-37]. 

3.  Electrochemistry  and  related  phenomena.  Side  reactions  inside  the 
battery  and  double-layer  capacitance  are  not  considered,  which 
is,  again,  a  common  assumption  [14,38,39]. 


Table  2 

PCM,  electrolyte,  and  outer  can  properties. 


Parameter 

Unit 

PCM  [1] 

Electrolyte  [22] 

Outer  can 

cp 

J(kgK)-1 

1980 

2055 

475 

k 

W(mK)-1 

16.6 

0.6 

44.5 

P 

kg  ITT3 

866 

1129.95 

7850 

L 

Jkg-1 

00 

X 

o 

- 

- 

Ts 

K 

325 

- 

- 

Ti 

K 

328 

- 

- 

$ 

- 

0.9  (graphite  alone) 

- 

0.8 
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Table  3 

Constants  in  Urefj,  dUrefti/dT,  and  cr;  expressions  [10,15]. 


Constant 

Unit 

Value 

Constant 

Unit 

Value 

Constant 

Unit 

Value 

m 

mV  K1 

344.1347 

P3 

mVK-1 

-26.0645 

- 

-3.0 

n2 

- 

-32.9633 

P4 

mVK-1 

12.7660 

V 

10 

n3 

- 

8.3167 

P5 

mVK-1 

4.3127 

01s 

- 

-2000 

n4 

- 

1 

P6 

- 

0.5715 

Vi 

v 

4.1983 

n5 

- 

749.0756 

P7 

mVK-1 

-0.1842 

¥2 

v 

0.0565 

ri6 

- 

-34.7909 

Ps 

- 

-0.5169 

^3 

- 

-14.5546 

n7 

- 

8.8871 

P9 

- 

0.0462 

^4 

- 

8.6094 

ns 

mVK-1 

-0.8520 

PlO 

mVK”1 

1.2816 

^5 

v 

-0.0275 

ng 

mVK-1 

0.3622 

Pll 

- 

-4.9916 

^6 

- 

0.9984 

nio 

mVK”1 

0.2698 

Pl2 

mVK-1 

-0.0904 

^7 

- 

0.4924 

a0 

S  m1 

1.0793  x  10“2 

Pl3 

- 

-20.9669 

^8 

- 

-1.9011 

d\ 

S  m2  mol-1 

6.7461  x  10“4 

Pl4 

- 

-12.5788 

^9 

v 

-0.1571 

a2 

S  m5  mol-2 

-5.2245  x  10-7 

Pl5 

mVK-1 

0.0313 

^PlO 

- 

-0.0474 

a3 

S  m8  mol-3 

1.3605  xlO-10 

Pl6 

- 

31.7663 

Vll 

v 

0.8102 

a4 

S  m11  mol-4 

-1.1724  x  10-14 

Pl7 

- 

-22.4295 

^Pl2 

- 

-40 

Pi 

mV  K"1 

-4.1453 

Oil 

V 

-0.16 

^Pl3 

- 

-0.1339 

P2 

mV  K1 

8.1471 

ot2 

V 

1.32 

- 

4.  Material  properties.  Uniform  distribution  of  active  materials  of 
uniform  size  in  the  electrodes  and  PCM  is  assumed.  Further, 
phase  change  for  the  PCM  does  not  occur  at  a  single  temperature 
but  rather  over  a  given  melting  range,  which  is  assumed  to  be 
captured  reasonably  well  with  a  linear  relation  between  latent 
heat  in  the  ‘mushy’  region  -  where  a  liquid  and  solid  phases 
coexist  -  of  the  PCM.  The  emissivity  of  PCM  is  taken  to  be  the 
emissivity  of  pure  graphite  alone. 

5.  Microscale  properties.  The  active  material  is  assumed  to  be  spher¬ 
ical;  i.e.,  we  only  need  to  consider  the  radial  direction  at  the 
microscale. 

6.  Contact  resistance.  The  contact  resistance  between  the  PCM  and 
the  battery  as  well  as  between  functional  layers  is  assumed  to 
be  negligible  as  a  first  approximation. 

The  governing  equations,  boundary  conditions  and  constitu¬ 
tive  relations  are  summarized  in  Appendices  A-D.  The  various 
parameters  and  constants  are  summarized  in  Tables  1-3. 


3.  Numerics 

The  commercial  finite-element  solver,  COMSOL  Multiphysics 
4.1  [40],  was  employed  to  solve  the  2D  model  after  importing  the 
geometry  from  AutoCAD  2011  [41  ],  in  which  the  2D  spiral-wound 
battery  with/without  PCM  was  created  with  the  helix  command 
by  setting  the  height  of  the  helix  to  zero  (z-direction).  The  wounds 
were  drawn  one  after  the  other  starting  from  the  center  (x  =  0,y  =  0) 
in  order  to  avoid  interference  of  the  various  layers  that  occurred 
when  all  the  wounds  were  drawn  in  a  single  step.  The  imported 
geometry  appeared  as  a  curve  in  COMSOL,  whence  it  had  to  be 
coerced  to  a  solid,  followed  by  splitting  the  formed  solid  into  sep¬ 
arate  subdomains  representing  the  various  functional  layers.  In 
addition,  the  narrowing  gap  of  liquid  electrolyte  between  the  jelly 
roll  and  the  outer  can,  as  illustrated  in  Fig.  If,  was  dealt  with  by 
treating  the  electrolyte  present  in  the  vicinity  of  the  region  where 
the  jelly  roll  touches  the  outer  can  as  part  of  the  latter. 

Linear  elements  were  implemented  for  all  dependent  variables: 
0S,  0/,  Ci,  cssurf ,  T  and  c™8',  the  direct  solver  UMFPACK  was  chosen 
as  linear  solver  with  a  relative  convergence  tolerance  of  10-3;  and 
solutions  for  all  models  were  tested  for  mesh  independence. 

Charge  and  discharge  currents,  iapp,  were  applied  at  the  respec¬ 
tive  boundaries  with  a  smoothed  Heaviside  function. 

Furthermore,  in  order  to  avoid  numerical  instabilities  due  to 
negative  values  in  the  current  density  when  the  local  state  of  charge 
(SOC)  approaches  unity  (i.e.,  fully  charged),  a  relational  operator 


Fig.  2.  (a)  Cell  voltage  vs  time  and  (b)  power  vs  time  for  galvanostatic  discharge  at 
various  C-rates  with  (symbols)  and  without  PCM  (lines). 
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was  introduced  for  the  exchange  current  density  and  local  SOC,  to 
ensure  that  the  former  does  not  become  negative  or  exactly  zero: 


here,  i0  is  the  exchange  current  density,  F  is  the  Faraday’s  constant, 
ko  is  the  reaction  rate  constant,  q  is  the  concentration  of  Li-ions 
in  the  electrolyte,  cssur^  is  the  concentration  of  Li-ions  on  the  sur¬ 
face  of  the  active  material,  c™ax  is  the  maximum  concentration  of 
Li-ions  in  the  active  material,  0pe  is  the  local  SOC  of  the  positive 
electrode,  and  c*  =  10-12  mol  m-3  is  a  negligible  concentration  cho¬ 
sen  such  that  it  does  not  affect  the  solution  at  leading  order;  i.e., 


t  =1 800  s 


t  =  3600  s 


Fig.  3.  Local  distribution  of  the  following  dependent  variables  at  t=  1 800  s  and  t  =  3600  s  during  discharge  at  1  C-rate:  SOC  of  positive  electrode  (a,  b),  SOC  of  negative  electrode 
(c,  d),  lithium  ion  concentration  in  the  electrolyte  (e,  f),  and  liquid  phase  potential  (g,  h). 
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The  relational  operators  are  introduced  only  for  the  positive 
electrode  since  the  latter  determines  the  capacity  of  the  considered 
battery;  more  on  this  later  in  the  discussion. 

All  computations  were  carried  out  on  a  workstation  with  two 
quad-core  processors  (3.2  GHz,  with  a  total  of  eight  processor  cores) 
and  a  total  of  64  GB  random  access  memory  (RAM). 


4.  Results  and  discussion 

In  the  following,  we  will  first  study  the  behavior  during  dis¬ 
charge  of  a  spiral-wound  lithium-ion  battery  without  a  passive 
thermal-management  system,  such  that  only  natural  convection 
and  radiation  provide  heat  transfer  with  the  surrounding.  (We  focus 
on  the  discharge  process,  since  charging  exhibits  a  similar  behav¬ 
ior.)  Thereafter,  we  compare  the  behavior  of  the  cell  with  a  PCM 
layer  surrounding  the  battery. 


41.  Discharge  and  power  curves 

We  start  with  the  global  behavior  in  terms  of  cell  voltage  and 
power  during  discharge  for  three  different  C-rates  of  the  battery, 
as  depicted  in  Fig.  2a  and  b.  Overall,  the  battery  discharge  occurs, 
as  expected,  first  gradually  with  decreasing  cell  voltage  and  power 
with  respect  to  time,  followed  by  a  sharp  drop  towards  the  end  of 
discharge.  From  the  numerical  point-of-view,  the  latter  is  reflected 
in  a  slowing  convergence  rate  and  requires  the  introduction  of  rela¬ 
tional  operators  as  introduced  in  Numerics  earlier.  Further,  we  note 
that  the  discharge  curves  are  more  or  less  the  same  for  a  battery 
with  and  without  PCM;  more  on  this  later  in  the  discussion. 


42.  Edge  and  geometry  effects 

In  a  spiral-wound  Li-ion  cell,  the  inner  and  outer  ends  of  the 
spiral-wound  jelly  roll  -  see  Fig.  1  d  and  e  respectively  -  are  exposed 
to  the  liquid  electrolyte  both  in  the  core  as  well  as  in  the  gap 
between  the  jelly  roll  and  the  outer  can.  Furthermore,  the  inner¬ 
most  layer  comprising  the  positive  electrode  is  only  in  contact  with 
the  electrolyte  in  the  core  and  a  current  collector  in  the  first  wound ; 
similar  for  the  negative  layer  in  the  outermost  wound,  which  is  only 
in  contact  with  a  current  collector  and  the  separator  layer.  These 
edge  and  geometrical  aspects  from  the  spiral-wound  geometry  can 
thus  be  expected  to  give  rise  to  localized  effects  and  deviations 
from  the  average  behavior.  This  is  indeed  the  case,  as  can  be  inferred 
from  Fig.  3,  for  the  SOC,  lithium-ion  concentration  in  the  electrolyte 
and  liquid-phase  potential  during  discharge  at  1  C-rate  halfway 
through,  t= 30  min,  and  at  the  end,  t  =  60  min,  for  a  lithium-ion  bat¬ 
tery  without  PCM.  Here,  two  main  features  are  apparent:  first,  is  the 
uniformity  in  the  interior  of  the  cell  and  second,  the  expected  devi¬ 
ations  that  occur  at  the  inner  and  outer  edges  of  the  spiral  wounds. 
For  example,  the  SOC  of  the  positive  electrode  at  the  inner  edge  of 
the  spiral  half-way  through  the  discharge  is  around  0.2  (Fig.  3a), 
which  is  close  to  the  initial  SOC  of  0.1 7 ;  at  the  outer  edge  of  the  spi¬ 
ral  and  in  the  interior,  the  SOC  is  around  0.9  and  0.55  respectively. 
There  is  thus  a  lack  of  lithium-ions  in  the  first  wound,  and  an  excess 
in  the  outer.  At  the  end  of  the  discharge,  the  SOC  approaches  1  for 
the  positive  electrode  except  at  the  inner  edge  of  the  spiral  where 
it  is  around  0.3  (Fig.  3b).  The  same  phenomena  are  observed  for  the 
negative  electrode  (albeit  reversed)  as  well  as  for  the  lithium-ion 
concentration  and  liquid-phase  potential.  These  findings,  in  turn, 
suggest  that  the  positive  electrode  at  the  inner  edge  of  the  spi¬ 
ral  and  the  negative  electrode  at  the  outer  edge  of  the  spiral  are 
not  completely  utilized  -  hence,  these  regions  can  be  left  uncoated 
while  manufacturing  a  battery. 


Time  /  s 

Fig.  4.  Time  history  of  (a)  the  average  battery  temperature  and  (b)  the  temperature 
difference  between  the  battery  core  and  the  surface  during  discharge  at  various  rates 
with  (dotted)  and  without  PCM  (continuous). 


4.3.  Heat  generation  and  thermal  behavior 

As  alluded  to  in  the  introduction,  thermal  management  is  key  in 
ensuring  not  only  a  safe  operation  but  also  to  improve  the  cycle-life 
of  a  lithium-ion  battery.  In  this  case,  the  battery  without  a  PCM  is 
only  cooled  through  natural  convection  and  radiation,  which  gives 
rise  to  a  significant  temperature  increase  inside  the  cell  during  dis¬ 
charge,  as  shown  in  Fig.  4a.  The  temperature  increases  by  around 
8  K,  20  K,  and  50  K  above  the  ambient  temperature  for  discharge 
rates  of  1  C,  2  C,  and  5  C,  respectively,  with  a  maximum  tempera¬ 
ture  of  350 1<  reached  at  the  end  of  discharge  at  the  highest  C-rate 
considered  here,  which  indicates  the  need  for  a  thermal  manage¬ 
ment  system.  The  radiation  from  the  cell  accounts  for  around  60% 
of  the  total  heat  transfer  to  the  ambient  as  compared  to  40%  by 
natural  convection.  The  temperature  differential,  AT,  between  the 
core  and  the  outer  can  of  the  battery  (Fig.  4b)  is  negligible  com¬ 
pared  to  the  overall,  average  temperature  inside  the  battery;  i.e. 
AT/Tavg~  10-3  <  1,  which  indicates  that  the  energy  transport  is 
mainly  limited  by  the  energy  exchange  with  the  ambient. 

The  rise  in  temperature  originates  from  the  electrochemical 
reaction  and  ohmic  heating  as  charge  is  passed  through  the  battery. 
The  total  heat  generation  rate  and  the  contribution  from  each  of  the 
heat  sources  (see  Eq.  (D.6))  during  discharge  at  1  C-rate  and  5  C-rate 
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are  shown  in  Fig.  5a  and  b  respectively.  At  a  1  C-rate,  the  reversible 
heat  generation  is  higher  than  the  other  sources,  amounting  for 
nearly  50%  of  the  total  heat  generation,  followed  by  ohmic  and  irre¬ 
versible  heating.  The  reason  can  be  traced  back  to  the  functional 
form  of  dUref  ildT  given  by  Eqs.  (D.7)  and  (D.8),  which  determines 
the  behavior  of  the  reversible  heating  with  time:  at  a  1  C-rate,  the 
reversible  heat  is  negative  initially  and  then  it  changes  sign  as  the 
discharge  proceeds,  whereas  the  other  heat  sources  are  always 
positive.  The  ohmic  heat  generation  starts  rising  during  the  ini¬ 
tial  period  of  discharge  due  to  the  currents  that  are  passed  through 
the  battery  and  then  almost  becomes  a  constant  for  the  rest  of  the 
discharge. 

At  a  5  C-rate,  however,  the  ohmic  heat  generation  is  the  highest 
contributor,  amounting  for  nearly  50%  of  the  total  heat  genera¬ 
tion,  followed  by  the  irreversible  and  the  reversible  heat  sources. 
The  ohmic  heat  generation  becomes  increasingly  larger  due  to  the 
increasing  current  flow  through  the  battery.  However,  a  drop  in 
the  ohmic  heat  generation  is  observed  after  the  initial  increase 
unlike  1  C  rate  where  it  remains  constant.  The  increase  in  battery 
temperature  is  mirrored  by  a  decrease  in  the  electric  resistance, 
whence  ohmic  heating  decreases  during  the  discharge;  in  the 
herein  derived  model,  the  electric  resistance  of  the  battery  is  cou¬ 
pled  with  the  temperature-dependent  ionic  conductivity  and  mass 
diffusivity  in  the  solid  phase.  In  contrast,  the  irreversible  heating  is 


a 


almost  constant  with  time  with  a  slight  increase  towards  the  end  of 
discharge. 

A  comparison  of  the  total  heat  generation  (Fig.  5a  and  b)  with 
the  power  delivered  by  the  battery  (Fig.  2b)  suggests  that  nearly 
10%  of  the  battery  power  is  lost  as  heat  at  a  1  C-discharge  rate;  in 
comparison,  25%  is  lost  through  heat  for  the  5  C-rate. 

Besides  noting  the  relative  magnitude  of  the  various  heat 
sources  during  discharge,  it  is  also  of  interest  to  see  which  layer  con¬ 
tributes  the  most  to  the  heat  generation  (Fig.  5c  and  d):  the  negative 
electrode  generates  most  of  the  heat,  amounting  to  nearly  60%  of 
the  total  heat  both  at  1  C  and  5  C-rates,  followed  by  the  positive  and 
then  by  the  remaining  layers.  The  reason  for  this  behavior  can  be 
found  in  the  ohmic  and  reversible  heating,  which  are  higher  in  the 
negative  electrode  than  in  the  positive  electrode  due  to  lower  ionic 
conductivity  and  the  thicker,  less  porous  nature  of  the  negative 
electrode  -  the  former  is  mainly  dependent  on  the  material  prop¬ 
erties  and  the  latter  are  design-adjustable  parameters.  The  heat 
generation  in  the  current  collector,  electrolyte  and  separator,  on 
the  other  hand,  only  comprise  ohmic  heating,  whence  they  remain 
almost  constant  throughout  the  galvanostatic  discharge. 

Finally,  we  address  the  local  temperature  distribution,  as  illus¬ 
trated  in  Fig.  6,  half-way  and  at  the  end  of  discharge  for  a 
1  C-rate.  Compared  to  the  other  dependent  variables  discussed 
earlier,  the  temperature  does  not  exhibit  any  edge  effects; 

b 


Fig.  5.  Time  history  of  heat  generation  by  various  sources  (a,  b),  total  heat  generation  and  heat  generation  in  various  layers  (c,  d)  during  discharge  at  1  C  and  5  C-rates  with 
(dotted)  and  without  PCM  (continuous). 
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instead,  the  temperature  distribution  is  near-to  axially  symmetric 
(z-axis). 

4.4.  Passive  thermal  management 

Thus  far,  we  have  seen  that  the  limiting  heat  removal  from 
the  cell  by  natural  convection  and  radiation  is  responsible  for  the 
overall  temperature  increase  and  that  the  temperature  distribution 
inside  the  cell  is  axially  symmetric  with  an  overall  temperature  gra¬ 
dient  that  is  negligible  compared  to  the  average  temperature  in  the 
battery.  One  approach  to  manage  the  thermal  envelope  could  be 
in  the  form  of  forced  or  mixed  external  convection  with  air/liquid 
through  active  thermal  management;  another  approach,  and  the 
one  we  pursue  here,  is  through  passive  thermal  management  with 
a  PCM  as  an  additional  layer  surrounding  the  battery  cell.  In  essence, 
the  disadvantage  of  doing  so  is  that  the  additional  layer  effectively 
adds  additional  resistance  to  heat  removal  from  the  battery;  the 
advantage,  however,  is  that  doing  so  increases  the  thermal  capaci¬ 
tance. 

Returning  to  Fig.  2,  we  note  that  the  discharge  curves  and  power 
characteristics  remain  unaffected  at  leading  order  when  a  PCM 
layer  with  thickness  of  1  mm  is  added  to  the  cell  -  which,  at 
first  sight,  would  suggest  that  the  passive  thermal  management 
is  ineffective  and  unnecessary.  In  contrast,  however,  the  average 
temperature  (Fig.  4a)  is  lowered  as  compared  to  the  same  cell  with¬ 
out  the  PCM,  whereas  the  average  temperature  gradient  inside  the 
cell  (Fig.  4b)  is  slightly  higher.  The  former  can  be  explained  by  the 
increase  in  the  thermal  capacitance  of  the  system  (battery  +  PCM) 
and  the  latter  by  the  increased  resistance  to  energy  transfer  out 
of  the  system.  Further,  the  PCM  has  not  yet  reached  its  melting 
temperature  under  the  discharge  rates  of  1  C  and  2  C,  whence  the 
advantage  of  cooling  through  phase-change  is  not  realized  in  these 
two  cases.  Under  a  discharge  rate  of  5  C,  however,  the  heat  gener¬ 
ated  by  the  battery  increases  its  temperature  to  the  melting  range 
of  the  PCM,  such  that  the  battery  is  around  18  K  cooler  at  the  end 
of  discharge  as  compared  to  without  PCM  -  a  substantial  decrease. 
Furthermore,  at  a  5  C-rate,  the  temperature  differential  with  PCM 
reaches  a  maximum  of  1.2  K  and  then  decreases  rapidly  owing  to 


the  drop  in  heat  generation  from  the  battery  and  the  loss  of  heat  to 
the  ambient  through  convection  and  radiation;  it  increases  again 
as  heat  generation  starts  to  rise. 

The  layer  with  PCM  also  affects  the  heat  generated  by  the 
battery  at  the  5  C-rate  whereas  at  a  1  C-rate,  there  is  no  notice¬ 
able  difference  in  heat  generation,  as  shown  in  Fig.  5.  Overall, 
the  total  heat  generation  increases  by  nearly  7%  with  PCM  as 
compared  to  without  at  a  5  C-rate.  This  corresponds  to  the  fact 
that  the  presence  of  PCM  keeps  the  battery  temperature  lower  as 
compared  to  the  battery  without  PCM,  which  in  turn  lowers  the 
ionic  conductivity  and  mass  diffusivity,  thus  increasing  the  resis¬ 
tance  of  the  battery  and  finally  leading  to  an  increase  in  ohmic 
heating. 

The  PCM  also  maintains  the  temperature  uniformity  within  the 
cell  unlike  forced-convection  cooling  [1,2],  for  which  appreciable 
temperature  differences  between  the  core  and  the  outer  surface  of 
the  battery  can  be  established.  The  temperature  distribution  inside 
the  battery  with  the  PCM  is  shown  in  Fig.  6c  and  d,  in  which  the 
temperature  rise  from  t = 320  s  to  t = 640  s  is  1 6 1<  without  PCM  and 
4  K  with  PCM  at  a  5  C-rate.  The  maximum  temperature  difference 
between  the  core  of  the  battery  and  the  outer  can  is  0.8  K  without 
PCM  and  1  K  with  PCM  and  this  increase  in  gradient  is  due  to  the 
additional  resistance  for  heat  transfer  to  the  ambient  as  mentioned 
above. 

The  state  of  the  PCM  can  readily  be  identified  as  solid,  liquid  or 
mushy  based  on  the  temperature  of  the  system.  When  it  is  in  the 
mushy  region,  another  parameter  is  needed  to  exactly  define  its 
state:  the  liquid  fraction.  For  1  C  and  2C-rates,  since  the  PCM  has 
not  yet  reached  its  melting  temperature  range,  the  liquid  fraction 
is  zero  whereas  at  a  5  C  rate,  the  wax  present  inside  the  graphite 
matrix  starts  melting,  but  does  not  reach  a  completely  liquid  state 
at  the  end  of  discharge  as  the  liquid  fraction  is  0.7  (Fig.  7),  indicat¬ 
ing  that  it  lies  in  the  two  phase  or  the  mushy  region.  The  volume  of 
the  PCM  in  the  layer  surrounding  the  battery  should  thus  be  cho¬ 
sen  so  that  it  is  not  completely  liquid  at  the  end  of  discharge  to 
ensure  cooling  throughout  the  discharge.  The  drop  at  around  360  s 
in  the  liquid  fraction  during  discharge  mirrors  the  drop  in  the  heat 
generation  inside  the  battery,  as  shown  earlier  in  Fig.  5b. 
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Fig.  7.  Liquid  fraction  of  PCM  vs  time  during  discharge  at  5  C-rate. 


With  this  model,  one  can,  e.g.,  optimize  the  thickness  of  the  PCM 
while  accounting  for  the  operating  conditions  of  the  battery. 

5.  Conclusions 


Appendix  A. 

A 1.  Governing  equations  (macroscale) 


The  mathematical  model  comprises  conservation  of  species  and 
charge  in  the  solution-  and  the  solid-phase  as  well  as  the  conser¬ 
vation  of  energy  [2,10,15,43]: 


V-i s  =  -J (pe,  ne,  cc), 

(A.1) 

V  •  i(  =  J(pe,  ne,  sp,  el), 

(A.2) 

s>  ft  +  V ' Nf  =  F  (pe>  ne’ sp’ el)’ 

(A3) 

dT 

pCp  +  V  q  =  () (pe,  ne,  sp,  cc,  el), 

(A.4) 

3  u 

P~r — b  V  •  q  =  O(PCM), 
ot 

(A.5) 

where  the  fluxes  are  given  by 

II 

i 

Q 

< 

(A.6) 

2  RTrreff 

i,  =  -ofVth  +  —p~  (1  -  t°)  V(ln  q). 

(A.7) 

z 

II 

1 

< 

n 

+ 

**  |  +  o 

(A.8) 

A  two-dimensional  coupled  thermal-electrochemical  model  for 
a  commercially  available  spiral-wound  Li-ion  battery  has  been  pre¬ 
sented  and  analyzed  for  discharge  in  terms  of  geometry  and  edge 
effects  as  well  as  in  terms  of  passive  thermal  management  with  a 
PCM.  The  reduction  in  dimensionality  was  justified  through  scal¬ 
ing  arguments  and  negligible  heat  flux  in  the  axial  direction  of  the 
battery. 

In  summary,  the  active  material  is  depleted  to  a  larger  extent 
at  the  outer  end  of  the  spiral  for  the  positive  electrode  and  at 
the  inner  end  of  the  spiral  for  the  negative  electrode  than  at  the 
remaining  parts,  where  it  is  utilized  uniformly.  Further,  reversible 
heat  generation  is  the  highest  contributor  at  lower  discharge  rates 
around  1  C-rate,  whereas  ohmic  heating  is  the  highest  contributor 
among  the  heat  generation  sources  for  higher  discharge  rates.  The 
findings  suggest  that  the  battery  design  can  be  optimized  in  order 
to,  for  example,  reduce  the  ohmic  heating  by  improving  the  elec¬ 
trolyte  conductivity  and  the  design-adjustable  parameters  like  the 
thickness  and  porosity  of  the  electrodes. 

A  passive  thermal  management  system  has  been  evaluated  by 
wrapping  a  PCM  around  the  battery,  which  reduces  the  average 
temperature  of  the  battery  for  higher  discharge  rates:  here,  around 
5  C-rate. 

The  model  can  be  extended  to  include  the  effect  of  current  col¬ 
lecting  tabs,  other  active  materials,  and  the  axial  dimension  if  the 
heat  flux  is  not  negligible  by  extruding  the  geometry.  Furthermore, 
thermal  runaway  can  be  simulated  by  including  the  additional 
reaction  heats  and  temperature  dependence  of  the  reaction  rates. 
Finally,  the  numerical  procedure  to  implement  the  geometry  and 
mathematical  formulation  could  be  automated  [42]  to  allow  for 
wide-ranging  parameter  studies  as  well  as  multi-objective  opti¬ 
mization  of  a  Li-ion  battery  cell  with  passive  thermal  management. 
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q  =  -kVT.  (A.9) 

In  the  above  equations,  is  is  the  solid-phase  current  density  J  is  the 
transfer  current  per  unit  volume,  crf^  is  the  effective  conductivity 
of  the  solid  phase,  0S  is  the  solid  phase  potential,  i/  is  the  solution 
phase  current  density,  is  the  effective  electric  conductivity  of 
the  solution  phase,  0/  is  the  solution  phase  potential,  R  is  the  uni¬ 
versal  gas  constant;  T  is  the  absolute  temperature,  F  is  the  Faraday’s 
constant,  is  the  transference  number  of  the  cation,  q  is  the  vol¬ 
ume  fraction  of  electrolyte  in  the  electrodes  and  the  separator,  and 
N i  is  the  molar  flux  of  the  cations;  p  is  the  effective  density  of  the 
various  functional  layers,  Cp  is  the  effective  specific  heat  capacity, 
q  is  the  conductive  heat  flux,  H  is  the  enthalpy  of  the  PCM,  Q  is 
the  heat  generation  per  unit  volume,  De^  is  the  effective  diffusion 
coefficient  in  the  liquid/solution  phase,  and  k  is  the  effective  ther¬ 
mal  conductivity;  the  electrolyte  is  a  binary  electrolyte  with  the 
concentration  of  the  electrolyte  q  defined  as  [30]  q  =  c+/v+  =  c_/v_, 
where  c+  and  c_  are  the  concentrations  of  the  cations  and  anions 
respectively,  and  v+  and  v_  represent  the  number  of  cations  and 
anions  produced  by  the  dissociation  of  one  mole  of  electrolyte.  (The 
governing  equations  are  applicable  in  the  layers  mentioned  inside 
parentheses.) 


Appendix  B. 


BA.  Governing  equations  (microscale) 


At  the  microscale,  lithium  diffuses  into  the  spherical  active 
material  particles  (Fig.  lg),  whence  the  conservation  of  lithium 
inside  the  active  material  can  be  stated  as 


3 cs  1  9  /  2n  9cs\ 

3 t=^dr\  Ds^FJ’ 


(B.l) 


where  cs{x,y,  r,  t)  is  the  concentration  of  protons  in  the  active  mate¬ 
rial  particle  of  the  electrode ;  t  represents  the  time,  Ds  is  the  diffusion 
coefficient  of  lithium  in  the  active  material,  and  r  is  the  radial  coor¬ 
dinate  inside  an  agglomerate.  Flere  we  employ  the  diffusion  length 
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approach  [44]  and  a  polynomial  approximation  [45],  such  that  the 
governing  equations  in  the  microscale  are  reduced  to 


^(cfr/-cD  = 

Is 

dcr_  3  ifm 

dt  F93  * 


ll 

F9 


(B.2) 

(B.3) 


cssurf  and  c“vg  are  the  volume-averaged  surface  and  average  con¬ 
centrations  of  Li-ion  in  the  active  material,  if  is  the  charge  transfer 
current  density  and  is  the  radius  of  the  active  material  in  the 
electrodes. 


Appendix  C. 

C.l.  Boundary  and  initial  conditions 

At  the  interface  I  between  the  electrode/electrolyte  or  the 
separator/electrolyte  (see  Fig.  Id  and  e  for  placement  of  roman 
numerals),  we  specify  continuity  for  the  energy  flux  as  well  as 
the  ionic  flux  of  lithium  ions  (ionic  current),  whereas  insulation 
is  specified  for  the  solid-phase  current: 

n •  is  =  0,  n-i/|4  =n-i/|L,  n-N/|4  =  n-N/|i_  (I),  (C.l) 

The  current  density  is  prescribed  at  the  positive  current  collector 
in  the  inner  end  of  the  spiral: 

n  •  is  =  -iapp  (II).  (C.2) 

At  the  current  collector/electrode  interfaces,  continuity  of  energy 
flux  and  solid-phase  current  is  specified;  insulation  is  specified  for 
the  ionic  flux  and  current: 

n  is|„4  =  n-islm  ,  n  q|„4  =  n  q|m_, 

nij  =  nN/  =  0  (III).  (C.3) 

At  the  electrode/separator  interfaces,  we  define  continuity  of 
energy  flux  and  ionic  flux  as  well  as  ionic  current  and  since  there  is 
no  flow  of  electrons  across  the  interface,  insulation  for  solid  phase 
current  is  defined. 

n  •  is  =  0,  n-i/|iV+  =n-ii\w_, 

n  q|IV+  =  n-q|IV_n-N/|IV+  «=  n-Nj|iv_  (IV).  (C.4) 


The  negative  current  collector  at  the  outer  end  of  the  spiral  is 
grounded: 

0s  =  O(V).  (C.5) 

At  the  electrolyte/can  interface,  there  is  continuity  of  energy  flux 
and  no  flow  of  ions: 

n  •  q|v4  =  n  •  q|Vi_ ,  n  i,  =  n  N,  =  0(VI).  (C.6) 

At  the  outer  surface  of  the  can,  both  convection  and  radiation  is 
considered: 

n  ■  q  =  t,(T  —  Ta)  +  tf{T4  -  T4)  (VII).  (C.7) 

When  the  battery  is  covered  with  PCM,  there  is  continuity  of  energy 
flux  from  the  can  to  the  PCM  and  on  the  outer  surface  of  the 
PCM  (VIII),  Newton’s  law  of  cooling  is  specified  along  with  radia¬ 
tive  energy  transfer  as  given  by  Eq.  (C.7).  In  Eqs.  (C.1)-(C.7),  n 
denotes  the  unit  normal  vector  for  a  given  boundary  or  interface, 
iapp  is  the  applied  current  density,  §  is  the  emissivity,  and  is 
the  Stefan-Boltzmann  constant.  The  battery  is  discharged  under 
galvanostatic  conditions  at  various  current  densities. 

At  t  =  0,  we  prescribe 

cs  =  c*vg  =  cl  c*  =  Cj°  (C.8) 


0s  = 


0  (ne) 
4>s  (pe) 


(C.9) 


4>i  =  4>°  (ne,  pe.  sp)  (C.10) 

T  =  T0  (C.ll) 


Appendix  D. 

DA.  Constitutive  relations  and  parameters 


The  source  term,  J,  is  expressed  as 

I  =  /  Aif  (ne’Pe)  mi) 

J  \0  (sp,  cc,  el)  ’  1  j 

where  A  is  the  specific  surface  area  for  the  faradaic  reaction  per 
unit  volume;  the  local  charge  transfer  current  density  is  given  by 
the  Butler-Volmer  equation  for  electrode  kinetics 


(D.2) 


Here,  i0  is  the  exchange  current  density,  aa  and  ac  are  the  anodic 
and  cathodic  transfer  coefficients  and  p  is  the  overpotential.  The 
specific  interfacial  area  is  related  to  the  particle  radius,  volume  void 
fraction  of  filler  €f,  polymer  matrix  ep,  and,  solution  phase  [10]  as 


A  = 


3(1  -€,-€f 

SH 


fp) 


(D.3) 


The  exchange  current  density  is  given  by 


i0  =  Fk0  (D.4) 

/<o  is  the  reaction  rate  constant  and  c™ax  is  the  maximum  lithium 
concentration  in  the  electrodes.  The  overpotential  is  defined  as 

=  uff  i,  i  =  ne,  pe  (D.5) 

in  which  Uef  is  the  effective  open  circuit  potential  of  the  elec- 
trode  with  respect  to  the  solid  lithium  electrode.  The  internal  heat 
generation  rate  Q.  [15]  is  defined  as 

Q=jri  +  af(y<t>s)2  +  ofiS/fr)2 


+  ' 


2  RTo\ 


f 


(1  -  t+)V(ln  q)  •  V0/,  i  =  ne,  pe 


(D.6) 


The  entropic  change  as  a  function  of  state  of  charge  for  the  negative 
electrode  is  expressed  as 

da 


ref,  ne  _  tM  exp(ri2^ne  +  ^3)  a  \  Q2  \ 

0= —  - ; - — Z - ; - \  +  n8^ne  +  ttgtfne  +  nio, 

dT  n4  +  n5  exp(n6t^ne  +  ^7) 
and  for  the  positive  electrode  as 


dU, 


ref,  pe 


3  T 


—  Pi  +  P2^pe  +  P3@pe  +  P4^pe  +  P5  exP(Pe3pe) 


+  p7  exp  - 


Ope  +  Ps 
P9 


+  P10  sin(pn6>Pe) 


(D.7) 


+  P12  sin(pis3pe  +  P14)  +  p  1 5  sin(p163Pe  +  P17),  (D.8) 


in  which  in  and  ip  [15]  are  constants  obtained  from  curve  fitting 
with  experimental  data;  3ne  and  3pe  denote  the  state  of  charge  of 
the  negative  and  positive  electrodes  respectively,  and  are  defined 
as 

surf 

One  =  0pe  =  “Sax  * 

Ls 


(D.9) 
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The  first  term  on  the  right-hand-side  (RHS)  in  the  heat  generation, 
Eq.  (D.6),  captures  the  irreversible  heat  generation  arising  due  to 
the  reaction  that  is  responsible  for  the  deviation  of  the  potential 
from  the  equilibrium  potential;  the  second  term  is  the  reversible 
heat  due  to  the  changes  in  entropy;  the  third  term  quantifies  the 
ohmic  heating  in  the  solid  phase;  and  the  last  two  terms  reflect 
the  ohmic  heating  in  the  solution  phase.  The  overall  state  of  charge 
of  the  battery  here  refers  to  the  local  state  of  charge  of  the  posi¬ 
tive  electrode  as  the  battery  capacity  is  limited  due  to  the  positive 
electrode.  The  electrode  capacity  is  defined  as 

q  =  w,(l  -e,-€f-  €p)PjCch,  i  =  ne,  pe,  (D.10) 

where  w,  represents  the  thickness  of  the  electrodes  and  Cth  is  the 
theoretical  capacity  of  the  electrode  material  (372  mAhg-1  for  the 
negative  electrode,  and  1 48  mAh  g-1  for  the  positive  electrode).  The 
total  battery  capacity  will  be  the  minimum  of  the  two  electrode 
capacities  and  in  our  case,  the  positive  electrode  has  the  minimum 
capacity. 

The  open-circuit  potential  of  the  two  electrodes  are  taken  from 
Doyle  et  al.  [10]:  for  the  positive  electrode,  we  write 

Uref, pe  =  ^Pl  +  ^2  tanh(<J330pe  +  ^4)  +  ^5  _  q — py  +  ^8^ 

+  %  exP  (^lO^pe)  +^Pll  exP(^Pl2(^pe  +  ^13))’  (D.ll) 
and  for  the  negative  electrode 

Uref, ne  =  ^1  +  ^2  exp(9t30ne)  +  ^4  exp(9t56>ne),  (D.12) 

where  iVl  and  Ep are  constants  obtained  by  curve  fitting  with  exper¬ 
imental  data. 

The  effective  (i.e.,  temperature-dependent)  open-circuit  poten¬ 
tial  of  an  electrode  is  approximated  by  a  first  order  Taylor-series 
expansion  around  a  reference  temperature  Tref  [38]: 

Uref,i  =  Uref,i  +  (j  ~  Tref )  *  =  ne,  pe  (D.13) 

The  effective  conductivities  and  diffusivity  are  defined  as 
erf  =  crse,  (D.14) 

af  =  a,eY,  (D.15) 

Df  =  Dley,  (D.16) 

where  y  is  the  Bruggemann  constant. 

The  electrolyte  conductivity  is  expressed  as  a  function  of  con¬ 
centration  c/  [10]: 


The  initial  values  of  the  solid  phase  and  the  liquid  phase  poten¬ 
tial  are  given  by 

<t>°s  =  Uref, pe(0°pe )  -  Ure/, ne(C ),  (D.20) 

0°  =  -lW,  ne(0  (D.21) 

where  0pe  and  denote  the  initial  states  of  charge  (SOC)  of 
the  positive  and  negative  electrodes  respectively.  The  solid-phase 
potential  difference  between  the  positive  current  collector  at  the 
inner  end  of  the  spiral  and  the  negative  current  collector  at  the 
outer  end  of  the  spiral  is  defined  as  the  cell  voltage,  Ecen\ 

Ecell  =  0S III  -  I V-  (D.22) 


The  current  density  at  which  the  battery  becomes  completely 
discharged  in  1  h  is  taken  as  the  1  C-rate.  (Complete  discharge  here 
implies  the  local  SOC  in  the  positive  electrode  reaches  one.) 

The  phase  change  for  the  material  considered  here  takes  place 
over  a  given  melting  range,  which  is  captured  the  following  func¬ 
tional  form,  H(T)  [43]: 


CPT 


H(T)  =  { 


CpT  T 


T-Ts 


-P  '  t,-Ts 
CpT  +  L 


for  T  <  Ts  (solid  region), 
for  Ts  <T  <  T/(mushy  region), 
for  T  >Ti  (liquid  region), 


(D.23) 


where  L  is  the  latent  heat  of  the  phase  change  and  Ts  and  7)  are  the 
start  and  end  temperatures  of  the  phase  change  respectively. 

The  average  temperature  for  the  battery  is  defined  as 


(D.24) 


where  Rb  is  the  radius  of  the  battery  and  Q  is  the  region  of  inte¬ 
gration  representing  the  entire  battery  as  shown  in  Fig.  lb.  The 
temperature  differential  is  defined  as 


AT  =  T |r=0  - 


1 

27 TRb 


(D.25) 


where  t  is  the  radial  coordinate  and  the  second  term  on  the  right 
hand  side  represents  the  average  temperature  at  the  outer  surface 
of  the  can.  The  total  heat  generation  and  the  heat  generation  in 
various  layers  is  defined  as 


<«>-'■// 


QdA 


(D.26) 


(f)  =  h  [  [ 

J  Jn 


QidA , 


i  =  ne,  pe,  sp,  el,  cc 


(D.27) 


4 

ct' = 

1=0 


(D.17) 


where  a{  are  constants. 

The  physical  properties  -  diffusion  coefficients  (both  solid  and 
liquid)  and  electrical  conductivity  (liquid)  -  are  dependent  on  tem¬ 
perature,  the  function  of  which  is  typically  written  in  the  form  of 
an  Arrhenius  expression  [15,38]: 


©(T)  =  ®(Tre/)  exp 


(D.18) 


where  @(T)  is  a  placeholder  for  a  temperature-dependent  property, 
Trej  is  a  reference  temperature,  and  Ea  @  is  the  activation  energy. 

The  diffusion  length,  /s,  for  the  spherical  electrode  particles  is 
estimated  as  [44] 

ls  =  j.  (D.19) 


where  h  is  the  height  of  the  battery  considering  only  the  jelly  roll 
without  the  top  cap  assembly. 
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